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This article makes the key observation that when using cyhndrical al- 
gebraic decomposition (CAD) to solve a problem with respect to a set of 
polynomials, it is not always the signs of those polynomials that are of 

^~H paramount importance but rather the truth values of certain quantifier 

^ free formulae involving them. This motivates our definition of a Truth 

CQ Table Invariant CAD (TTICAD). We generalise the theory of equational 

constraints to design an algorithm which will efficiently construct a TTI- 
CAD for a wide class of problems, producing stronger results than when 
using equational constraints alone. The algorithm is implemented fully in 

\l Maple and we present promising results from experimentation. 
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1 Introduction 

Cylindrical algebraic decompositions (CADs) are a key tool in real algebraic ge- 
ometry, both for their original motivation, solving quantifier eUmination prob- 
lems, but also for use in many other applications ranging from robot motion 
planning [22' etc.] to programming with complex functions [12] etc.]. Tradi- 
tionally CADs are produced sign-invariant to a given set of polynomials, (the 
signs of the polynomials do not vary on the cells of the decomposition) . How- 
ever, this gives far more information than required for most problems. The idea 
of a truth invariant CAD (the truth of a formula does not vary on each cell) was 
defined in [2] for use in simplifying CADs. The key contribution of this paper 
is an approach to construct CADs which are truth invariant without having to 
first build a sign-invariant CAD. Actually, we directly build CADs which are 
truth table invariant, (the truth values of various quantifier free formulae do not 
vary). 



We present an algorithm to efficiently produce TTICADs for a wide class 
of problems, utilising the theory of equational constraints ^9j. The algorithm 
goes further than equational constraints by allowing the creation of smaller 
CADs in a wider variety of cases; for example disjunctive normal form where 
each individual conjunction has an equational constraint but no single explicit 
equational constraint is present for the formula. The problem of decomposing 
complex space according to a set of branch cuts for the purpose of algebraic 
simplification ([2TJ etc.]) is of this case. 

1.1 Background on CAD 

We briefly remind the reader about the theory of CAD, first proposed by Collins 
in [g. 

Definition 1. A Tarski formula -F(xi, . . . , x„) is a Boolean combination (A, V, ^) 
of statements about the signs, (— 0, > 0, < 0, but therefore 7^ 0, > 0, < as 
well), of certain integral polynomials fi{xi, . . . ,a;„). We use QFF to denote a 
quantifier free Tarski formula. 

CAD was developed as a tool for the problem of quantifier elimination over 
the reals: given a quantified Tarski formula 

Qk+lXk+l ■ ■ ■ Q„XnF{xi, ...,x„) (1) 

(where Qi £ {V, 3} and F is a QFF), produce an equivalent QFF V'(a;i, ■ . . , Xk)- 
Collins proposed to decompose M" cylindrically such that each cell was sign- 
invariant for all fi occurring in F. Then tp would be the disjunction of the 
defining formulae of those cells c^ in K'^ such that (fTl) was true over the whole 
of Cj, which is the same as saying that (fTl) is true at any one "sample point" of 

Ci- 

Collins' algorithm has two phases. The first, projection, applies a projection 
operator repeatedly to a set of polynomials, each time producing another set 
in one fewer variables. Together these sets contain the projection polynomials. 
These are then used in the second phase, lifting, to build the CAD incrementally. 
First M is decomposed into cells which are points and intervals corresponding 
to the real roots of the univariate polynomials. Then M^ is decomposed by 
repeating the process over each cell using the bivariate polynomials at a sample 
point. The output for each cell consists of sections (where a polynomial vanishes) 
and sectors (the regions between). Together these form a stack over the cell, 
and taking the union of these stacks gives the CAD of M^ . This is repeated until 
a CAD of M" is produced. 

To conclude that a CAD produced in this way is sign-invariant we need 
delineability. A polynomial is delineable in a cell if the portion of its zero set 
in the cell consists of disjoint sections. A set of polynomials are delineable in 
a cell if each is delineable and the sections of different polynomials in the cell 
are either identical or disjoint. The projection operator used must be defined 
so that over each cell of a sign-invariant CAD for projection polynomials in r 
variables, the polynomials in r -I- 1 variables are delineable. 



The output of a CAD algorithm depends on the variable ordering. We usu- 
ally work with polynomials in Z[a;i, . . . ,a;„] with the variables, x, in ascending 
order (so we first project with respect to Xn and continue to reach univariate 
polynomials in a;i). The main variable of a polynomial (mvar) is the greatest 
variable present with respect to the ordering. 

Major directions of work since 1975 includes the following: 

1. Improvements in Collins' main algorithms by |17l and many others]. These 
have focussed on reducing the projection sets required as discussed further 
later. 

2. Complexity theory of CAD 0(13]. 

3. Partial CAD, introduced in [TT], where the structure of F is used to lift 
only when required to deduce ip. 



4. The theory of equational constraints, [THl HOI IS discussed in Section 2.1 
This is related to the previous direction but differs by using more efficient 
projections. 

5. CAD via Triangular Decomposition '8]: a radically different approach for 
computing a sign-invariant CAD which is used for Maple's inbuilt CAD 
command. 

1.2 TTICAD 

We define a new type of CAD, the topic of this paper. 

Definition 2. Let <& — {(f>i}l^i be a list of QFFs. We say a cylindrical alge- 
braic decomposition V is a Truth Table Invariant CAD for $ (TTICAD) if the 
Boolean value of each 4>i is constant (either true or false) on each cell ofT>. 

A full sign-invariant CAD for the set of polynomials occurring in the formulae 
of $ would clearly be a TTICAD. However, we aim to produce an algorithm 
that will construct smaller TTICADs for certain $. We will achieve this using 
the theory of equational constraints (first suggested in [TU] with the key theory 
developed in |il9j). 

Definition 3. Suppose some quantified formula is given: 

</>* = {Qk+lXk+l) ■ ■ ■ {QnXn)<l>{^)- 

where the Qi are quantifiers and (j) is quantifier free. An equation f — is 
called an equational constraint of (p* if f ~ is logically implied by (j> (the 
quantifier-free part of (j)* ). Such a constraint may be either explicit or implicit. 

We suppose that we are given a formula list <& in which every QFF (/),j has a 
designated explicit equational constraint ft — 0. We will construct TTICADs by 
generalising McCallum's reduced projection operator for equational constraints 
(as in [19]) so that we may make use of the equational constraints. 



Figure 1: The polynomials from Section 1.3 



1.3 Worked Example 

We will provide details for the following worked example. 
Consider the polynomials: 

/i := a;^ + y^ - 1 .gi := xy - ^ 

/2 := {x - 4)2 + (y - 1)2 - 1 52 := (x - 4)(y - 1) - i 

which are plotted in Figure [T] We wish to solve the following problem: find the 
regions of M^ where the formula 

$ :- (A = A gi < 0) V (/2 = A 32 < 0) 

is true. Assume that we are using the variable ordering y )~ x (so the 1- 
dimensional CAD is with respect to x). 

Both Qepcad [3] and Maple 16 [8] produce a full sign-invariant CAD for 
the polynomials with 317 cells. At first glance it seems that the theory of 
equational constraints [111 [IHl IS] is not applicable here as neither /i = nor 
/2 = is logically implied by $. However, while there is no explicit equational 
constraint we can observe that /1/2 = is an implicit constraint of $. Using 
Qepcad with this declared gives a CAD with 249 cells. Later, in Section [273| 
we demonstrate how a TTICAD with 105 cells can be produced. 



2 Projection Operators 

2.1 Equational Constraints 

We use two key theorems from McCallum's work on projection and equational 
constraints. Both theorems use CADs which are not just sign-invariant but 
have the stronger property of order-invariance. A CAD is order-invariant with 
respect to a set of polynomials if each polynomial has constant order of vanishing 
within each cell. 

Let P be the McCallum projection operator [17], which produces coefficients, 
discriminant and cross resultants from a set of polynomials. We assume the 
usual trivial simplifications such as removal of constants, exclusion of entries 
identical to a previous entry (up to constant multiple), and using only the 



Figure 2: Graphical representation of Theorem |2] 




necessary coefficients. Recall that a set ^4 C Z[x] is an irreducible basis if the 
elements of A are of positive degree in the main variable, irreducible and pairwise 
relatively prime. The main theorem underlying P follows. 

Theorem 1 (jl8)). Let A be an irreducible basis in Z[x] and let S be a connected 
submanifold o/M"^^. Suppose each element of P{A) is order-invariant in S. 
Then each element of A either vanishes identically on S or is analytic delineable 
on S, (a slight variant on traditional delineability, see il8^ )- The sections of 
A not identically vanishing are pairwise disjoint, and each element of A not 
identically vanishing is order-invariant in such sections. 

The main mathematical result underlying the reduction of P in the presence 
of an equational constraint / is as follows. 

Theorem 2 ( 19 ). Let /(x),(7(x) be integral polynomials with positive degree 
in Xn, let r{xi, . . . ,Xn-i) be their resultant, and suppose r ^ 0. Let S be a 
connected subset ofM"~^ such that f is delineable on S and r is order-invariant 
in S. Then g is sign-invariant in every section of f over S. 

Figure [2] gives a graphical representation of the question answered by Theo- 
rem [2] Here we consider polynomials f{x, y, z) and g{x, y, z) of positive degree 
in z whose resultant r is non-zero, and a connected subset 5 C M^ in which r 
is order-invariant. We further suppose that / is delineable on S (noting that 
Theorem 1 with n — "i and A — {/} provides sufficient conditions for this). We 
ask whether g is sign-invariant in the sections of / over S. Theorem [2] answers 
this question affirmatively: the real variety oi g either aligns with a given section 
of / exactly (as for the bottom section of / in Figure [2]), or has no intersection 
with such a section (as for the top). The situation at the middle section of / 
cannot happen. Theorem [2] thus suggests a reduction of the projection operator 
P relative to an equational constraint / = for the first projection step, as in 

m. 



2.2 A Projection Operator for TTICAD 

In [19] the central concept is that of the reduced projection of a set A of integral 
polynomials relative to a nonempty subset E oi A and it is an extension of this 
which is central here. For simplicity in [TU], the concept is first defined for the 
case when A is an irreducible basis and by analogy we start with a similar special 
case. Let A — {Ai}l^^ be a list of irreducible bases Ai and let £ = {Ei}^^-^ be a 
list of nonempty subsets Ei C- Ai. Put A — [j^^i Ai and E = (Jj^^ Ei (we will 
use the convention of uppercase Roman letters for sets and calligraphic letters 
for sequences). 

Definition 4. We define the reduced projection of A with respect to £ , denoted 
by P£{A), as follows: 



^£(-4):=ULi^s.(A0URcs><(f) (2) 



whe 



PeAA) - PiE^) U {res,„(/,ff) \feE,,gG A,,g^ E,} 
Res^ {£) = {res,„ (/, /) | 3z, j : / e E, J & Ej,t < jj ^ /} 

In Scction [3.1| we build Algorithmlllto apply the reduced projection operator 
for less special input sets by considering contents and irreducible factors of 
positive degree. 

Definition 5. The excluded projection polynomials of {Ai,Ei) are those in 
P{A) but excluded from P£{A): 

ExcIPb, (A,) := P(AO \ Pe, (A) (3) 

= {coeSs{g),disc^^{g),Tes^^{g, g) \ g,g e A,\Ei,g ^ g}. 

The total set of excluded polynomials, denoted ExclP£(A), consists of all the 
ExclP^;. (Ai), along with the cross resultants of gi with all of Aj for i ^ j- 

The following theorem is an analogue of Theorem 2.3 of _19J, and provides the 
foundation for our algorithm in Section [3. 1[ 



Tiieorem 3. Let S be a connected submanifold of'W^~^. Suppose each element 
of P£{A) is order invariant in S. Then each f € E either vanishes identically 
on S or is analytically delineable on S, the sections over S of the f £ E which 
do not vanish identically are pairwise disjoint, and each element f £ E which 
does not vanish identically is order-invariant in such sections. 

Moreover, for each i, with 1 < i < t, every g £ Ai\ Ei is sign-invariant in 
each section over S of every f £ Ei which does not vanish identically. 

Proof. The crucial observation is that P{E) C P£{A). To see this, recall equa- 
tion (pi) and note that we can write 

P{E)=[j^PiE,)URes''{£). 



We can therefore apply Theorem [T] to the set E and obtain the first three 
conclusions immediately. 

There remains the final conclusion to prove. Let i be in the range 1 < i < t, 
let g £ Ai\ Ei and let f £ Ei; suppose / does not vanish identically on S. 
Now veSx^ifi g) € P£{A)., and so is order-invariant in S by hypothesis. Further, 
we already concluded that / is delineable. Therefore by Theorem [2j g is sign- 
invariant in each section of / over S. D 

In the following section we can use Theorem [3] as the key tool for our imple- 
mentation of TTICAD, so long as the equational constraint / does not vanish 
identically on the lower dimensional manifold, S. When working with a poly- 
nomial / considered in r variables that vanishes identically at a point a e M'"^^ 
we say that / is nullified at a. 

Remark 4. It is clear that the reduced projection P£{A) will lead to fewer (or 
the same) projection polynomials than the full projection P. One may con- 
sider instead using the reduced projection Pe{A) of 119], (with E = UiEi and 



A — UiAi as above). In the context of Section 1.2 this corresponds to using 
W^ fi as an implicit equational constraint for a single formula. Note that Pg {A) 
also contains fewer polynomials than Pe{A) in general since Pe{A) contains all 
resultants res{f,g) where f € Ei,g G Aj (and g ^ E), while P£{A) contains 
only those with i = j (and g ^ Ei). 

2.3 Worked Example 

In Section |3] we will discuss how to use these results to define an algorithm 
for TTICAD. First we illustrate the potential savings with our worked example 
from Section [T31 

In the notation introduced above we have: 

Ai := {/i,.gi}, E, := {A}; A^ := 1/2,52}, E2 := {A}. 
We construct the reduced projection sets for each (j)i, 
PB,(^i) = {x2-l,a;4„a.2_,__i.}^ 

Pe2 (^2) = {x^ ~8x + 15, x^ - 16x^ -I- 95x^ - 248x + ^ } 
and the cross-resultant set 

Res'' {£) = {resy(/i, A)} = {GSx^ - 272a; -I- 285}. 

PsiA) is then the union of these three sets. In Figure Is] we plot the poly- 
nomials (solid curves) and identify the 12 real solutions of PsiA) (solid vertical 
lines). We can see the solutions align with the asymptotes of the /s and the 
important intersections (those of /i with gi and /2 with 52)- 

If we were to instead use a projection operator based on an implicit equa- 
tional constraint /1/2 = then in the notation above we would construct Pe{A) 



Figure 3: The polynomials from the worked example along with the solutions 
to the projection sets. 




Figure 4: Magnified region of Figure |3] 




from A — {/i, /2, 51,52} and E — {/i,/2}- This set provides an extra 4 solu- 
tions (the dashed vertical lines) which align with the intersections of /i with 
g2 and /2 with gi. Finally, if we were to consider P{A) then we gain another 
4 solutions (the dotted vertical lines) which align with the intersections of 51 
and 52 and the asymptotes of the gs. In Figure [4] we magnify a region to show 
explicitly that the point of intersection between /i and 51 is identified in P£{A)., 
whereas the intersection points of 52 with both /i and 51 are ignored. 

Hence the 1-dimensional CAD produced using PeiA) has 25 cells compared 
to 33 when using Pe{A) and 41 when using P{A). However, it is important 
to note that this reduction is amplified after lifting (using Theorem a\ and and 
Algorithm [1]). The full dimensional TTICAD has 105 cells, the CAD invariant 
with respect to the implicit equational constraint has 249 cells and the full 
sign-invariant CAD has 317. 



3 Implementation 

3.1 Algorithm Description and Proof 

We describe carefully Algorithm [l] This will create a TTICAD of E" for a 
list of QFFs, $ = {4'i}i=ii in variables x = xi ^ X2 ^ ■ ■ ■ ^ Xn where each 
cj)i has a designated equational constraint fi = of positive degree. We use a 
subalgorithm CADW, fully specified and validated in [18 . The input of CADW is: 



r, a positive integer and A, a set of r-variate integral polynomials. The output 
is a Boolean w which if true is accompanied by an order-invariant CAD for A 
(a list of indices / and sample points S) . 

Let Ai be the set of all polynomials occurring in 0^, put Ei = {fi}, and let 
A and £ be the lists of the Ai and Ei, respectively. Our algorithm effectively 
defines the reduced projection of A with respect to £ using the special case 
of this definition from the previous section. The definition amounts to using 
*P := C U Pj^{B) for P£{A), where C is the set of contents of all the elements 
of all the Ai, B is the list {Bi}*=i, such that Bi is the finest squarefree basis 
for the set prim(^i) of primitive parts of elements of Ai which have positive 
degree, and J" is the list {Fi}l^j^, such that Fi is the finest squarefree basis for 
piim{Ei). (The reader will notice that this notation and the definition of PsiA) 
is analogous to the work in Section 5 of [l9].) 

We shall prove that, provided A and £ are well-oriented as in Definition [6) 
the output of Algorithm [T] is indeed a TTICAD for $. Note that this condition 
is specialised and new, introduced for this paper. Its requirement is due to both 
the use of CADW from fl8| and the introduction of our new reduced projection 
operator. 

We first recall the more general notion of well-orientedness from [18] . A set 
A of n-variate polynomials is said to be well oriented if whenever n > 1, every 
/ £ prim(A) is nullified by at most a finite number of points in M"^^, and 
(recursively) P{A) is well-oriented. The Boolean output of CADW is false if the 
input set was not well-oriented in this sense. Now we define our new notion of 
well-orientedness for the set lists A and £ defined above, and hence <&. 

Definition 6. We say A is well oriented with respect to £ (and that $ is well 
oriented j if whenever n > 1, every constraint polynomial fi is nullified by at 
most a finite number of points in M"~^, and PsiA) (hence *P in the algorithm) 
is well-oriented in the sense of \18\j . 

Tiieorem 5. The output of Algorithm[^is as specified. 

Proof. We must show that when $ is well-oriented the output is a Truth Table 
Invariant CAD, (each (j)i has constant truth value in each cell of D), and FAIL 
otherwise. 

If the input was univariate then it is trivially well-oriented. The algorithm 
will construct a CAD V of M^ using the roots of the irreducible factors of the 
constraint polynomials (steps [5] to [6|. At each 0-cell all the polynomials in each 
(j)i trivially have constant signs, and hence every (pi has constant truth value. In 
each 1-cell no constraint polynomial has a root, so every (pi has constant truth 
value false. 

Now suppose n > 1. If *p is not well-oriented in the sense of [H] then CADW 
returns w' as false. In this case the input $ is not well oriented in the sense 
of Definition |6] and Algorithm [l] correctly returns FAIL. Otherwise, *P is well- 



oriented and at step 13 we have w' = true. Further, /' and S" specify a CAD, 
V , order-invariant with respect to *p. Let c, a submanifold of K""^, be a cell 
of v. 



Algorithm 1: TTICAD Algorithm 



Input : A list of quantifier-free formulae $ = {4>i}i=i in variables 

xi, . . . ,Xn- Each (/)j has a designated equational constraint 

/. = 0. 
Output: Either • V : A TTICAD of M" for $ (described by lists / and S 

of cell indices and sample points, respectively); or 

• FAIL: If $ is not well oriented (Def. |6]). 

for i = 1 . . . i do 
I Set Ei <— {fi}. Compute the finest squarefree basis Fi for prim(_Ei); 

Set F i^ Ul^j^F,; 
if n = 1 then 

Isolate in (/, S) the real roots of the product of the polynomials in F; 
return / and S for D; 
else 

for i = I . . .t do 

Extract the set Ai of polynomials in </>, ; 

Compute the set Ci of contents of the elements of Af, Compute 

the set Bi, the finest squarefree basis for prim(j4j); 

Set C ^ uU^a, B ^ (B,)Li and J" ^ (i^,)Li ; 
Construct the projection set: *P •(— C U Pj^{B) ; 
Attempt to construct a lower-dimensional CAD: 
u/,/',S" ^CADW(n- 1,^3); 
if w' = false then 
I return FAIL (*:p not well oriented); 

7^0; 5^0; 

for each cell c E V do 

Lc ^ {}; 

for i = 1 , . . . t do 

if fi is nullified on c then 
if dim(c) > then 
I return FAIL ($ not well oriented); 
else 
[_ L^^ L^UB,; 

else 

^ Lc^LcU F,; 

Lift over c using Lc'- construct cell indices and sample points for 
the stack over c of the polynomials in Lc, adding them to / and S] 

return / and S for V; 



1 

2 

3 
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9 
10 

11 
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13 
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Suppose first that the dimension of c is positive. If any constraint polynomial 
fi vanishes identically on c then $ is not well oriented in the sense of Definition [6] 
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and the algorithm correctly returns FAIL at step 22 Otherwise, we know that 
$ is certainly well-oriented. Since no constraint polynomial fi vanishes then no 
element of the basis F vanishes identically on c cither. Hence, by Theorem [3) 
applied with A = B and £ = T , each element of F is delineable on c, and the 
sections over c of the elements of F are pairwise disjoint. Thus the sections and 
sectors over c of the elements of F comprise a stack S over c. Furthermore, 
Theorem p^ assures us that, for each i, every element of Bi \ Fi is sign-invariant 
in each section over c of every element of Fi. 

Let 1 < i < t. Consider first a section a of the stack E. We shall show 
that (^i has constant truth value in a. Now the constraint polynomial ji is a 
product of its content cont(/i) and some elements of the basis Fi. But cont(/i), 
an element of *P, is sign-invariant in the whole cylinder c x M which includes a. 
Moreover all of the elements of Fi are sign-invariant in cr, as noted previously. 
Therefore ji is sign-invariant in cr. If /^ is positive or negative in a then 0^ has 
constant truth value false in a. 

Suppose that fi = throughout a. It follows that a must be a section of some 
element of the basis Fi. Let g € Ai \ Ei he a non-constraint polynomial in Ai. 
Now, by the definition of Bi, we see g can be written as g = cont{g)h^^ ■ ■ ■ h^'' 
where hj G Bi,pj e N. But cont(g), in *p, is sign-invariant in the whole cylinder 
c X K including a. Moreover each hj is sign-invariant in a, as noted previously. 
Hence g is sign-invariant in a. (Note that in the case where g does not have 
main variable Xn then g = cont(g) and the conclusion still holds). Since g 
was an arbitrary element of Ai\ Ei, it follows that all polynomials in Ai are 
sign-invariant in a, and hence that (j)i has constant truth value in a. 

Next consider a sector a of the stack E, and notice that at least one such 
sector exists. As observed above, cont(/i) is sign-invariant in c, and fi does not 
vanish identically on c. Hence cont(/i) is non-zero throughout c. Moreover each 
element of the basis Fi is delineable on c. Hence the constraint polynomial fi 
is nullified by no point of c. It follows from this that the algorithm does not 
return FAIL during the lifting phase. It follows also that fi^O throughout a. 
Therefore (jji has constant truth value false in a. 

It remains to consider the case in which the dimension of c is 0. In this case 
the roots of the polynomials in the lifting set Lc constructed by the algorithm 
determine a stack S over c. Each (jji trivially has constant truth value in each 
section (0-cell) of this stack, and the same can routinely be shown for each sector 
(1-cell) of this stack. D 

Remark 6. When the input to Algorithm^ is a single QFF then it produces a 
CAD which is invariant with respect to the sole equational constraint. This may 
be shown using the results of ]19j alone. However, we note that AlgorithmUjis 
actually more efficient in the lifting stage than the modified QEPCAD algorithm 
discussed in ]19^ since the lifting set excludes some non- equational constraint 
input polynomials. 

Algorithm [l] and Definition [6] have been kept conceptually simple to aid 
readability. However in practice the algorithm may sometimes be unnecessarily 
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cautious. In [3], several cases where non-well oriented input can still lead to an 
order-invariant CAD are discussed. Similarly here, we can sometimes allow the 
nullification of an equational constraint on a positive dimensional cell. 

Lemma 7. Let fi be an equational constraint which vanishes identically on a 
cell c E D' constructed during Algorithmln If all polynomials in ExclPe- (A^) 
are constant on c then any g £ Ai\ Ei will be delineable over c. 

Proof. Suppose first that Ai and E^ satisfy the simplifying conditions from Sec- 
2.2 Rearranging (Is]) we see 



tion 



P{A,) = PEAA^) U Exc1Pb,(A, 



However, given the conditions of the lemma, this is equivalent (after the 
removal of constants which do not affect CAD construction) to Pg. (yl^) on c. 
So here P{Ai) is a subset of PsiA) and we can conclude by Theorem fl] that all 
elements of Ai vanish identically on c or are delineable over c. 

In the more general case we can still draw the same conclusion because 
P{A,) = C, U Pf, {B,) U ExclPi.. {B,) C <p. D 

Hence we can use Lemma [7] to safely extend step [24] to also apply in such 
cases. In particular, we can allow equational constraints fi which do not have 
main variable Xn in such cases. We have included this in our implementation 



discussed in Section 3.3 In theory, we may be able to go further and allow step 
24 to apply in cases where the polynomials in ExclP^;. (A^) are not necessarily 
all constant, but have no real roots within the cell c. However, identifying such 
cases would require answering a separate quantifier elimination question, which 
may not be trivial. 

3.2 TTICAD via the ResCAD Set 

In Algorithm fl] the lifting stage (steps 16 to 27 1 varies according to whether 
an equational constraint is nullified. When this does not occur there is an 
alternative implementation of TTICAD which would be simpler to introduce 
into existing CAD algorithms. Define the ResCAD Set of $ as 

nm = i?U ULi {res,„(/,3) \f£E,,gE A,,g i E,} . 

Theorem 8. Let A = {Ai)\^^ be a list of irreducible bases Ai and let £ = 
{Ei)l^^ be a list of non-empty subsets Ei C Ai. For the McCallum projection 
operator P , \17l we have: 

P{TZm)^Pe{A). 

The proof is straightforward and so omitted here. 

Corollary 9. // no fi is nullified by a point in M"^^ then inputting ??.($) into 
any algorithm which produces a sign- invariant CAD using McCallum's projec- 
tion operator, will result in the TTICAD for $ produced by Algorithm^ 



12 



Hence Corollary|9]gives us a simple way to compute TTICADs using existing 
CAD implementations, such as Qepcad, but this cannot be applied as widely 
as Algorithm [T] 

3.3 Implementation in Maple 

There are various implementations of CAD available but none guarantee order- 
invariance, required for proving the validity of our TTICAD algorithm. Hence 
we needed to construct our own implementation to obtain experimental results. 
We built an implementation of McCallum projection, so that we could repro- 
duce CADW and modified the existing stack generation commands in Maple 
from [8] so they could be used more widely. Together these allowed us to fully 
implement Algorithm [T] The CAD implementation grew to a Maple package 
ProjectionCAD which gathers together algorithms for producing CADs via pro- 
jection and lifting to complement the existing CAD commands in Maple which 
use triangular decomposition, giving the same representation of sample points 
using regular chains. For further details (along with free access to the code) see 

3.4 Formulating a Problem for TTICAD 

When formulating a problem for TTICAD there may be choices for the input, 
such as choosing which equational constraint to designate in a QFF when more 
than one is present. Other possibilities include choosing whether conjunctions 
of formulae should be split into separate QFFs. Usually it will be preferable 
to minimise the number of QFFs, but if for example a designated equational 
constraint has many intersections with another polynomial which could be ig- 
nored by using separate QFFs, then the cost of the extra polynomials in the 
projection set may be outweighed by the complexity of those removed. Hence 
it is worth taking care in how we formulate the TTICAD. A simple problem of 
the form 

A - A /2 = A ffi < A 52 < 

has six acceptable choices for the composition of <&. 

We have started exploring heuristics for choosing the best composition. The 
metric sotd (sum of total degrees) as defined in [14' may be used to approximate 
the complexity of polynomials. We first considered using sotd(*P) and found 
that while it was fairly well correlated with the number of cells produced by 
Algorithm IT] it was not always fine enough to separate compositions leading 
to TTICADs with significantly different numbers of cells. Hence we prefer a 
stronger heuristic, sotd(*p U P(*P))) where P is the complete set of projection 
polynomials obtained by repeatedly applying P. 

For the problems in Section |4] we used the QFFs imposed by the disjunc- 
tions of formulae using this heuristic to choose which equational constraints are 
designated when there was a choice. For these problems the heuristic computa- 
tion time was negligible compared to the overall time, but for larger problems 
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this would not be the case. Work on heuristics is ongoing with a more detailed 
report available in ^. 

4 Experimental Results 

4.1 Description of experiments 

Our timings were obtained on a Linux desktop (3.1GHz Intel processor, 8.0Gb 
total memory) with Maple 16 (command hne interface), Mathematica 9 
(graphical interface) and Qepcad-B 1.69. For each experiment we produce 
a CAD and give the time taken and number of cells (cell count). The first is an 
obvious metric while the second is crucial for applications performing operations 
on each cell. 

For Qepcad the options +N500000000 and +L200000 were provided, the 
initialization included in the timings and explicit equational constraints declared 
when present with the product of those from the individual QFFs declared 
otherwise. In Mathematica the output is not a CAD but a formula constructed 
from one 2A, with the actual CAD not available to the user. Cell counts for 
the algorithms were provided by the author of the Mathematica code. 



TTICADs are calculated using our implementation described in Section 3.3 
which is simple and not optimized. The results in this section are not presented 
to claim that our implementation is state of the art, but to demonstrate the 
power of the TTICAD theory over the the conventional theory, and how it can 
allow even a simple implementation to compete. Hence the cell counts are of 
most interest. 

The time is measured to the nearest tenth of a second, with a time out (T/0) 
set at 5000 seconds. When F occurs it indicates failure due to a theoretical rea- 
son such as not well-oriented (in either sense). The occurrence of Err indicates 
an error in an internal subroutine of Maple's RegularChains package, used by 
ProjectionCAD. This error is not theoretical but a bug, beyond our control. 

We considered examples originating from [7 . However these problems (and 
most others in the literature) involve conjunctions of conditions, chosen as such 
to make them amenable to existing technologies. These problems can be tack- 
led using TTICAD, but they do not demonstrate its full strength. Hence we 
introduced new examples, denoted with a f , which are adapted from [7] to have 
disjuncted QFFs. 

Two examples came from the application of branch cut analysis for simpli- 
fication. These problems require a decomposition according to branch cuts of 
the form f — A g < 0, and then go on to test the validity of a simplification 
on each cell, iSl, etc.]. We need to consider the disjunction of the branch cuts 
making such problems suitable for Algorithm 111 We included a key example 
from Kahan ^6] , along with the problem induced by considering the validity of 
the double angle formulae for arcsin. Finally we considered the worked example 
from Section [T73| and its generalisation to three dimensions. Note that A and B 



following the problem name indicate different variable orderings. Full details for 
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all examples can all be found in the CAD repository [23] available freely online 
at http : //opus . bath . ac . uk/29503. 

4.2 Results 

We present our results in Table [l] For each problem we give the name used in 
the repository, n the number of variables, d the maximum degree of polynomials 
involved and t the number of QFFs used for TTICAD. We then give the time 
taken and number of cells produced by each algorithm. 

We first compare our TTICAD implementation with the sign-invariant CAD 
generated using ProjectionCAD with McCallum's projection operator [15]. Since 
these use the same architecture the comparison makes clear the benefits of the 
TTICAD theory. The experiments confirm the fact that the cell count for TTI- 
CAD will always be less than or equal to that of a sign-invariant CAD produced 
using the same implementation. Ellipsef A is not well-oriented in the sense of 
[18] . and so both methods return FAIL. Solotarefff A and B are well-oriented 
in this sense but not in the stronger sense of Definition [6] and hence TTICAD 
fails while the full sign-invariant CADs can be produced. The only example 
with equal cell counts is CoUisionf A in which the non-equational constraints 
were so simple that the projection polynomials were unchanged. Examining the 
results for the worked example and its generalisation we start to see the true 
power of TTICAD. In 3D Example A we see a 759- fold reduction in time and a 
50- fold reduction in cell count. 

We next compare our implementation of TTICAD with the state of the 
art in CAD: Qepcad 3 , Maple f5] and Mathematica [23] HI]. Mathe- 
MATICA is the quickest, however TTICAD often produces fewer cells. We note 
that Mathematica's algorithm uses powerful heuristics and so actually used 
Grobner bases on the first two problems, causing the cell counts to be so low. 
When all implementations succeed TTICAD usually produces far fewer cells 
than Qepcad or Maple, especially impressive given Qepcad is producing 
partial CADs for the quantified problems, while TTICAD is only working with 
the polynomials involved. For CoUisionf A the TTICAD theory offers no benefit 
allowing the better optimized alternatives to have a lower cell count. 

Reasons for the TTICAD implementation struggling to compete on speed 
in general are that the Mathematica and Qepcad algorithms are largely im- 
plemented directly in C, have had far more optimization, and in the case of 
Mathematica use validated numerics for lifting [5T. However, the strong 
performance in cell counts is very encouraging, both due its importance for ap- 
plications where CAD is part of a wider algorithm (such as branch cut analysis) 
and for the potential if TTICAD theory were implemented elsewhere. 

5 Conclusions 

We have defined Truth Table Invariant CADs, which can be more closely aligned 
to the needs of problems than traditional sign-invariant CADs. Theorem [3| ex- 
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tended the theory of equational constraints allowing us to develop Algorithm [T] 
to construct TTICADs efBciently for a large range of problems. The algorithm 
has been implemented in Maple giving promising experimental results. TTI- 
CADs in general have less cells than full sign-invariant CADs using the same 
implementation and we showed that this allows even a simple implementation 
of TTICAD to compete with the state of the art CAD implementations. It is 
anticipated that future implementations of TTICAD could be far better opti- 
mized leading to lower times for the same cell counts. We also note that the 
benefits of TTICAD increase with the number of QFFs in a problem and so 
larger problems may be susceptible to TTICAD when other approaches fail. 

We hope that these results inspire other implementations of TTICAD, with 
Corollary [9] showing a particularly easy way to adapt existing CAD implemen- 
tations. 

5.1 Future Work 

There is scope for optimizing the algorithm and extending it to allow less restric- 
tive input. Lemma[7|gives one extension that is included in our implementation 
while other possibilities include removing some of the caution implied by well- 
orientedness, analogous to \A^. Also, work developing heuristics for composing 
the input is underway in [T]. 

Of course, the implementation of TTICAD used here could be improved in 
many ways, but perhaps more desirable would be for TTICAD to be incorpo- 
rated into existing state of the art CAD implementations. In particular, we 
would like to use the existing Maple CAD commands [8 but this requires first 
understanding when they give order-invariance, a key question currently under 
consideration. We see several possibilities for the theoretical development of 
TTICAD: 

• Can we apply the theory recursively instead of only at the top level? For 
example by widening the projection operator to conclude order-invariance, 
as in gO]- 

• Can we apply TTICAD to forms of QFF other than "one equality and 
other items"? For example, can we generalise the theory of bi-equational 
constraints? 

• Can we make use of the ideas behind partial CAD to avoid unnecessary 
lifting once the truth value of a QFF on a cell is determined? 

• Can anything be done when $ is not well oriented? 

• Can we implement the lifting algorithm in parallel? 
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